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Abstract 

In the present paper we study subsolutions of the Dirac equation. The known examples of such 
subsolutions are the Weyl neutrino and the Majorana neutrino. It is shown that the Dirac equation 
for a massive non-interacting particle admits another decomposition into two separate equations. 
It is suggested that such particles can be identified with weakly interacting massive particles, other 
than neutrinos and of non-exotic nature, which can be candidates for dark matter. 
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I. INTRODUCTION 



The Dirac equation is a relativistic quantum mechanical wave equation formulated by 
Paul Dirac in 1928 providing a description of elementary spin-| particles, such as electrons 
and quarks, consistent with both the principles of quantum mechanics and the theory of 
special relativity [lj . The Dirac equation is 0, d, 0] : 

YPt$> = m*, (1) 

where m is the rest mass of the elementary particle, p M = i-^ are the components of 
the four-momentum operators defined in natural units (c = 1, H — 1). Tensor indices are 
denoted with Greek letters, jd = 0,1,2,3 and we shall always sum over repeated indices. 
The 7's are 4x4 anticommuting Dirac matrices: 7^7^ + 7^7^ = 2g flu I where g stands for 
the Minkowski space-time metric tensor: g^ v = diag (1, —1,-1, —1) and / is a unit matrix. 



The wave function is a bispinor, i.e. consists of two two-component spinors £, 77: ^ = 




In the present paper we study subsolutions of the Dirac equation. The known examples 
of such subsolutions are the Weyl neutrino and the Majorana neutrino, discussed in the 
next Section. In Section 3 the Dirac equation for a massive particle is decomposed into 
two separate equations. Lorentz covariance of these equations is demonstrated in Section 
4 where a symmetry of subsolutions is also found. Finally, possibility that such particles 
can be identified with weakly interacting massive particles, other than neutrinos and of 
non-exotic nature, is considered. 



II. WEYL AND MAJORANA NEUTRINOS AS SUBSOLUTIONS OF THE 
DIRAC EQUATION 



In the m = case it is possible to obtain two independent equations for spinors £, rj 

i7°7 1 7 2 7 3 



by application of projection operators Q± = | (1 ± 7 s ) to Eq.Q since 7 5 = - ; " — " ; 



anticommutes with 7 M p M : 

Q±'fp l & = 'fpAQ**) = °- (2) 
In the spinor representation of the Dirac matrices [3J] we have 7 s = diag (—1,-1,1,1) and 

thus Q-ty = ^ g j , = ^ |j j and separate equations for £, rj follow: 

(p° + ~a ■ -p) 77 = 0, (3a) 
{p° - ? ■ /) £ = 0, (3b) 

where a denotes the vector built of the Pauli matrices. Equations ([3]) are known as the 
Weyl equations and are used to describe massless left-handed and right-handed neutrinos. 
However, since the experimentally established phenomenon of neutrino oscillations requires 
non-zero neutrino masses, theory of massive neutrinos, which can be based on the Dirac equa- 
tion, is necessary 0,00|. Alternatively, a modification of the Dirac or Weyl equation, called 
the Majorana equation, is thought to apply to neutrinos. According to Majorana theory 
neutrino and antineutrino are identical and neutral Although the Majorana equations 
can be introduced without any reference to the Dirac equation they are subsolutions of the 
Dirac equation Q. 



2 



Indeed, demanding in ([I]) that \I/ = C^> where C is the charge conjugation operator, 
C^f = i7 2 \l>*, we obtain in the spinor representation £ = —ia 2 r]*, r\ = 2<T 2 £*and the Dirac 
equation ([!]) reduces to two separate Majorana equations for two-component spinors: 



It follows from the condition \? = that Majorana particle has zero charge built- 
in condition. The problem whether neutrinos are described by the Dirac equation or the 
Majorana equations is still open [1, @, @] . 

III. MASSIVE SUBSOLUTIONS OF THE FREE DIRAC EQUATION 
A. The spinor representation 

The free Dirac equation ([TJ) in the spinor representation of 7 matrices reads: 



(p° + ~~a ■ 'p') 77 = —ima 2 rj*, 
(p — a ■ p) £ = +ima 2 ^*. 



(4a) 
(4b) 



(p° + p 3 ) rji + (p 1 — ip 2 ) r\2 = m ? 1 1 
(p 1 + ip 2 ) % + (p° — p 3 ) r\2 = m £ 2 I 
(p° — p 3 ) £ x + (— p 1 + ip 2 ) £ 2 = mrji 1 
(— p 1 — ip 2 ) t; 1 + (p° + p 3 ) £ 2 = mr)2 J 



with \& = £ 2 , rji, r]^) T 0] (see also 0] for full exposition of spinor formalism). 
For m^Owe can define new quantities: 



(p° +p 3 ) Vi = m^fjj , (p 1 - ip 2 ) % = m([ 2) , 
(p 1 + ip 2 ) Vi = «^( 2 i) , (p° ~P 3 )V2 = m i 2 (2) , 



(6a) 
(6b) 



where we have: 




(7a) 
(7b) 



In spinor notation = ^J 1 , f£ 2) = > ffi) 



Equations ([5]) can be now written as: 



(p° + p 3 )r]i = m£L ' 
(p 1 - zp 2 ) 772 = mq 2) 
(p 1 + zp 2 ) r/i = rn&q 



(p° - p 3 ) r/ 2 = m£ ( 2 2) > . 



(p° - p 3 ) (efi) + ejjj) + (V + *p 2 ) fe } + £( 2 2) 1 = mr/i 




(9a) 
(9b) 
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Taking into account the identities 
the following two equations: 



we can finally write equations (jHJ) as a system of 



Due to the identities 



(p° + p 3 ) % 
(p 1 + ip 2 ) rji 
(P° ~ P 3 ) &) + (-p 1 + W 2 ) £ 



(P 1 - 
(p° 



(i) 

ip 2 ) rj 2 = 
P 3 ) V2 = 




i| equations ffTUl) . ffTTj) can be cast into form: 
(p° + p 3 ) rji = m£} x) 



m 

(p 1 + ip 2 ) r]i 



(p° ~ P 3 ) + (-P 1 + ip 2 ) t 



(i) 



(P 1 + IP 2 ) Cf; 



= mrji 




(i) 



(p 1 - ip 2 ) % 
(p° -p 3 )% 

{p° -p 3 )im + (-P l + w 2 )i 



mi 2 



(2) 



'(2) 

- zp 2 ) + (p° + p 3 ) £ ( 2 2) = m?7 2 



(2) 



(2) 





(10) 



(11) 



(12) 



(13) 



B. Arbitrary representation 

Let us define the following set of mutually commuting projection operators: 



Pi 

Pi 



1 (3- 7 5 - 7° 7 3 + z 7 V ) , 
1 (3- 7 5 + 7 °7 3 - *7V) , 



(14a) 
(14b) 
(14c) 
(14d) 



^3 = |(3+ 7 5 + 7 7 3 + nV), 
P 4 = I (3+ 7 5 - 7 °7 3 - ^7 1 7 2 ) , 

where 7 5 = — 27°7 1 7 2 7 3 , introduced by Corson 0] (note that in [9] operators q = 1 — Pi 
were defined and another convention for the 7 matrices, 7^7^ + 7^7^ = 2<5 A " y J, was used). In 
the spinor representation of 7 matrices we have Pi = diag (1, 1, 1, 0), P 2 = diag (1, 1, 0, 1). 

It is now easily checked that for m 7^ Eqs. (fT2]) . (TLB"]) can be written in the spinor 
representation of the Dirac matrices as 



7 M p M Pi%) = mFif(i), 

7^P2^(2) = ™P 2 #(2), 



(15a) 
(15b) 



respectively, where 



(i) 



£( 2 i) 

m 

\ V2 I 



(2) 



m 

V V2 ) 



Pl*(l) + P 2 ^( 2 ) = 



(16) 
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It is important that Eqs. (jl5ap . (115bj) are valid for arbitrary representation of the Dirac 



matrices since the projection operators (fl4j) are defined in arbitrary representation. 

Equations (flQj) . (TlTI) and identities (Q can be also given representation independent form. 
Indeed, from (j!5ap . (115bl) the equations 

P(fPrP& = mPi®, (i = 1, 2) (17a) 
(l-POyW^ . (i = 1, 2) (17b) 

which are equivalent to Eqs. (fID|) . (TTTT) and ([HD in the spinor representation, follow, respec- 
tively. Let us note that the identities (117bj) are the direct consequence of Eqs. f|T5l) for 



arbitrary representation of the Dirac matrices. 



IV. LORENTZ COVARIANCE AND SYMMETRIES OF SUBSOLUTIONS 

In what follows we shall use the approach of Bjorken and Drell [2j to study covariance 
of Eqs. (|T5|) . The postulate of Lorentz covariance means that the two equations below are 
equivalent: 

iT-^r ~ ™) %W W) = 0, (19a) 



where k — 1,2. Since all representations of 7 matrices are unitarily equivalent [10] we may 
write Eq. fll9ap as: 

^T^7-m)p^'(x') = 0. (19b) 



We start with transformation of the term P^ (x) in Eq. ( TTHTl : 

P k m (x) = PkS- 1 (a) (x') = S- 1 (a) [S (a) P^ 1 (a)] (x') = 

= S- 1 (a) (x') , (20) 

where S (a) is a similarity matrix depending on some set of parameters a and: 

P f h = S(a)P k S- 1 (a). (21) 

Substituting fl2DD into Eq.(JTSD we get: 



5 (a) (V^r - m) (a) P^' (*') = 0, (22) 
S (a) Ufa"^ - m\ S' 1 (a) Pp? {x>) = 0, (23) 
(iS (a) t^ 1 (a) a%^p - m\ (x') = 0. (24) 
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We now demand that: 

5 (a) y*^- 1 (a) a*; = 7". (25) 
Equations fl25|) have the following solutions: 

S = exp (-Ilog^P") , <V = 5 (7^ - 7^) , (26) 

where constant matrices define a Lorentz transformation 0]. 

Substituting Eq. (l25|) into (I21I) we obtain Eq. (119b|) and this completes the proof that 
Eqs.( TT8|) . (1221) are equivalent. It thus follows that Eqs. (|T5|) are covariant. 

There is a symmetry connecting equations (I15al) . (115bl) . To show this let us notice that 
it is possible to carry out Lorentz transformation to cast equations (I15al) . (115bl) into the 
following forms: 

(7°Po-7 1 Pi-m)P 1 % ) (a;) = 0, (27) 
(7% - 7V - m) P 2 ^ (2 ) (x) = 0, (28) 

respectively, since the Lorentz transformations in the x°x 3 and x 1 x 2 planes, are obtained 
from d2SD as S03 = e"^ 03 {I 30 = -I 03 = -1), S 12 = e^ ai2 {I 12 = -I 21 = -1), respectively 
[2|, commute with Pj's. Now it turns out that the unitary matrix V: 

V = ^7 2 7 3 , (29) 

commutes with 7 , 7 1 and transforms Pi <=^ P2 and thus transforms Eqs. f[2Tl) . (T2SJ) one into 
another in this special frame of reference. 



V. DISCUSSION 

It was demonstrated that, besides Weyl and Majorana neutrinos, there are also other 
subsolutions of the Dirac equation, c.f. Eqs. ffTol) . These new subsolutions are covariant and 
are valid for arbitrary representation of the Dirac matrices. 

Derivation of Eqs. (|T2|) . (|T3|) in the spinor representation of the Dirac matrices was 
based on the identities (Q. Analogous identities (117bj) are fulfilled in the case of arbitrary 
representation of the Dirac matrices. These identities are a consequence of commutation of 
components of four-momentum operator in the free case. It follows that in the interacting 
case, these equations cannot be derived from the corresponding Dirac equation. Moreover, 
the mass in Eqs. (1151) cannot be zero and these equations are not chiral. It thus follows that 
equations (1TB} describe some massive, spin-|, non-interacting and non-chiral particles. It 
seems natural to identify these particles with weakly interacting massive particles (WIMPs), 
other than neutrinos, which are candidates for dark matter, the existence of which was 
recently confirmed A strong candidate for WIMP dark matter is exotic particles, 



possibly the neutralino [121 . |13| |. On the other hand, theory based on Eqs. (115]) does not 
invoke exotic matter. 

Let us note finally that equations ( fT2l) . ( Tl3l) and their 3x3 analogues were also found in 



the Kemmer-Duffm-Petiau theory of spin and 1 mesons [14| . 
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